Site-selective nuclear magnetic relaxation time in a superconducting vortex state 
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The temperature and field dependences of the site-selective nuclear spin relaxation time Ti around 
vortices are studied comparatively both for s-wave and d-wave superconductors, based on the micro- 
scopic Bogoliubov-de Gennes theory. Reflecting low energy electronic excitations associated with the 
vortex core, the site selective temperature dependences deviate from those of the zero-field case, and 
Ti becomes faster with approaching the vortex core. In the core region, T" 1 has a new peak below 
the superconducting transition temperature T c . The field dependence of the overall Ti(T) behaviors 
for s-wave and d-wave superconductors is investigated and analyzed in terms of the local density 
of states. The NMR study by the resonance field dependence may be a new method to probe the 
spatial resolved vortex core structure in various conventional and unconventional superconductors. 

PACS numbers: 74.60.Ec, 74.25.Jb, 76.60.Pc 



I. INTRODUCTION 

Much attention has been focused on vortex mat- 
ter physics both of hiah T c cuprates and of conven- 
tional superconductors .Era There are several experimental 
methods to extract the vortex properties, ranging from 
thermodynamic measurements such as specific heat to 
various transport experiments such as thermalconduc- 
tivity under field. We can divide them into two cat- 
egories: One of the methods such as scanning tunnel- 
ing microscopy (STM) probes the electronic structure of 
the vortex, namely, low-lying excitations around a core. 
The other method probes the magnetic field distribution 
in the vortex state, examplified by small angle neutron 
diffraction, or muon spin resonance (/iSR) experiments. 
These two methods are complimentary to each other. We 
definitely need to increase the list of the method and fur- 
ther refine the theoretical analysis of experimental data 
to extract useful information about the vortex properties. 

Among various experimental methods mentioned 
above, the nuclear magnetic resonance (NMR) experi- 
ments is unique because it provides us two kinds of in- 
formation. Namely, the nuclear resonance spectrum re- 
flects the magnetic field distribution and the longitudi- 
nal relaxation time T\ probes electronic excitations in 
the vortex state through its T-dependence. The NMR 
experiment has been playing a vital role in distinguish- 
ing between s-wave and d-wave pairing symmetries in 
this respect, i.e. via T-dependence of T\. The power 
law Tj" 1 cx T 3 (T 5 ) behavior is taken as definitive ev- 
idence for a line (point) node in the gap structure of 
unconventional superconductors. This conclusion comes 
from a simple power counting of the density of states 
(DOS) at the Fermi level as a function of the energy 
E: N(E) cx E (E 3 ) for a line (point) node in a bulk 
superconductor at zero field. However, actual NMR ex- 
periments are performed under applied fields in a mixed 
state. The contributioiL.of the vortex core is inevitably 
included in their dataJj'B Usually, T\ is measured by se- 
lecting the resonance frequency at a most intensive sig- 



nal in the resonance spectrum. However, the resonance 
spectrum reflects information ofj-ijiternal magnetic field 
distribution of the vortex latticeErQ as mentioned above. 
Thus, by choosing the resonance field, we can specify the 
spatial position to detect the NMR signal. The signal at 
the maximum (minimum) cutoff comes from the vortex 
center (the furthest) site. The signal at the logarithmic 
singularity of the resonance field comes from the saddle 
points in the field distribution. By studying the position 
dependence of T\ around vortices through the resonance 
frequency dependence, we can clarify the detail of the 
vortex contribution in the NMR experiments. It helps us 
in the analysis of the standardized procedure extracting 
the gap symmetry. 

Low-lying excitation spectra around a vortex are not 
fully understood both experimentally and theoretically. 
The related problems are as follows. In the s-wave super- 
conductors, the effect of the quantized energy level will 
appear in the quasi-particle stateHrJ In the (i-wave case, 
the low energy state around the vortex core extends out- 
side the core due to the node of the superconducting gap, 
and shows, the viT-like DOS relation (H is an applied 
field) .E3~ll3 We also need to estimate the quasi-particle 
transfer between vortices (such as the path of the transfer 
and its amplitude) to study the dHvA-ascillation or trans- 
port phenomena in the mixed state.lHjO The excitation 
around the core plays a fundamental role in determining 
physical properties of superconductors. Several recent 
theories based on the microscopic Bogoliubov-de Gennes 
(BdG) equatiosi investigate electromagnetic properties of 
mixed stater.pftn jB^cmmection with the low-lying vortex 
excitations.BE3ll3'E2rE3 In high T c cuprates, the existence 
or non-existence of localized core excitations in <i-wave 
pairing case is actively debated. Theoretical study sug- 
gested the zero-energy peak in the d-wave case—instead 
of the quantized energy level in the s-wave case.Ejilj On 
the contrary, the STM experiments reported quantized 
energy level with large gap in YBCO (Ref. 23) and also 
in BSCCO (Ref. ^), and no peak within the supercon- 
ducting gap in BSCCO (Ref. E5h. A part of reasons of 
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the debate is due to limited experimental methods which 
directly probe the spatially resolved core structure. So 
far, the STM was only a method to detect it as mentioned 
before. A large number of thermodynamic or transport 
measurements probe spatially averaged quantities. 

Here we propose a novel spatially resolved means, that 
is, vortex imaging to see electronic excitations associated 
with a vortex core by using NMR, and demonstrate how 
the T-dependence of T\ is site-sensitive. Through this 
analysis, we are able to produce a spatial image of the 
low-lying excitation spectrum around a core. A similar 
idea of the NMR imaging is actually tested exjpftiimcn- 
tally in high T c materials-, by Slichter's groupErcZl and 
theoretically proposed,! 2 ^r. 9 ! .and also in spin-Peierls sys- 
tem CuGeC>3 by HorvaticEj 

In order to analyze the NMR data and propose suitable 
NMR experiments, we perform a model calculation of T± 
to demonstrate how careful NMR experiment is valuable. 
The other purpose of this paper is to warn a pitfall when 
obtaining the conclusion of the nodal gap structure from 
the T-dependence of Tj by performing the saddle point 
NMR. The position dependence of the NMR signal in the 
s-wave caseavas theoretically studied under some approx- 
imations ffl Here, we calculate it microscopically from 
the wave functions obtained by self-consistently solving 
the BdG equation for the extended Hubbard model in 
the s- and <i-wave cases. 

After giving the formulation of the problem in Sec. 
0, we show the results of basic vortex properties in this 
formulation in order to yield a coherent physical picture 
and to facilitate the later understanding of our results in 
Sec. IIL In next section the T-dependence of T\ is ana- 



lyzed both for s-wave and d-wave superconductors. The 
field-dependence of T\ is calculated in Sec. |v|. The final 
section is devoted to conclusion and discussions] 



II. FORMULATION 



A. Bogoliubov-de Gennes equation on lattice 



The BdG theory and its equation are one of 
the most fundamental frameworks in the theory of 
superconductivity.Ej In principle, the solution of BdG 
equation should give all static properties of a type II su- 
perconductor under an applied field, which we concern. 
Here we consider the BdG equation defined on a lattice 
instead of continuum space. The latter approach is use- 
ful when the order parameter is described by the s-wave 
symmetry, but difficult to treat non-local higher angular 
momentum states such as <i-wave. In contrast the lattice 
BdG theory is relatively easy to treat such a case and 
suitable for treating the s-wave and d-wave cases on an 
equal footing. 

In terms of the eigen-energy E a and the wave func- 
tions u a {vi), v a (vi) at i-site, the BdG equation for the 



extended Hubbard model defined on the two-dimensional 
square lattice is given by 



u a (v ) ( _ 



where 



with 



U,j = U,j exp [i — — / A(r) ■ dr 



(2) 



(3) 



(4) 



and the on-site interaction U, the chemical potential /i 
and the flux quantum <j)Q. The transfer integral U.j = t 
and the nearest neighbor (NN) interaction Vij = V for 
the NN site pair and Tj, and otherwise t% j = Vij = 0. 
The vector potential A(r) = x r in the symmetric 
gauge. The self-consistent condition for the pair poten- 
tial is 



Aij = -i^u Q (r J )<(r J )tanh( J B Q /2T). 



(5) 



This BdG equation is self-consistently solved by fallowing 
the numerical method of Wang and MacDonald.Ej 

In the following we examine comparatively the s-wave 
and d-wave symmetry cases. The s-wave pair potential 
is given by 



A s (r,) = UA U . 
The d x 2 _ J/ 2-wave pair potential is given by 
V 



(6) 



Ad(ri) = j(A £)I + A.^ - Ay :l - A_ u ) (7) 



with 



A±i ti = A M± ,sexp[i— 



(ri+r 4± e)/2 



A(r)-dr]. (8) 



We consider the square vortex lattice case where near- 
est neighbor vortex is located at the 45° direction from 
the a axis. This vortex lattice configuration is suggested 
for d-wave superconductors, or s-wave .supaxonductors 
with fourfold symmetric Fermi surface l 13 H 39 l'P°l The unit 
cell in our calculation is the square area of sites where 
two vortices are accommodated. Then, the magnetic field 
is determined as 



H 



2^0 
(cN r y 



(9) 



where c is the atomic lattice constant. For typical high 
T c cuprates, HN r=30 corresponds to an order of 15 Tesla. 
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We consider the area of N% unit cells. By introducing the 
quasi-momentum of the magnetic Bloch state, 

we set 

u a (r) = u Q (r)e ikr , v a (r) = £ Q (r)e ik - r . (11) 

We solve Eq. (0) within a unit cell under the given k. 
Then, a is labeled by k and the eigen- values obtained by 
this calculation within a unit cell. 



B. Boundary condition 

We impose the periodic boundary condition given by 
the symmetry for the translation 



R = rnui + nu.2 



(12) 



with m and n being integers, and Ui and 112 are unit 
vectors of the vortex lattice, i.e., 



g(x, x')=Tj2 e-^" (T - T,) 5(x, x', u n ), (19) 



the thermal Green functions with the Fermionic imagi- 
nary frequency ui n — 2nT(n + 5) are written as 



gu (x,x^„) = ^ — — 

a 

.912 (x,x ,U n ) = > — — 

z — ' IWn - E a 

a 

iui n — E a 



52l(x,x',LJ n ) = ^ 



(20) 
(21) 
(22) 



u a (r + R) = u a (r)e 



ix(r,R)/2 



(13) 



522 x,x ,w n ) = > — — . 23) 



The derivation of this form of the thermal Green func- 
tions is given in Appendix |a|. 



D. Local density of states 



v a (r + R) = v a (r)e 



-i X (r,R)/2 



(14) 



Here. 



27T 27T 

x(r, R) = - — A(R) • r - nmn + — (H x r ) ■ R (15) 

in the symmetric gauge when the vortex center is located 
at 



1 , 

r = r Q + -(ui + u 2 ). 



(16) 



The phase factor0 in Eq. (||) is needed to satisfy the 
translational relation 



A d (r + R) = A d (r)e i ^ R ). 



C. Thermal Green function 



(17) 



To understand the behavior of the position-dependent 
Ii(r), we also consider the local density of states (LDOS) 
at r. This is evaluated by using the thermal Green func- 
tions as 

NUE, r) = --Inu7n(r, r, iw -> E + i-q) (24) 

7T 

for the up-spin electron contributions, and 



N L (E,r) = -Img 22 (r,r,iu;-> E + irj) (25) 

7T 



for the down-spin electron contributions. Then, the 
LDOS is given by 

N(E,r)=Ni(E,r)+N l (E,r) 

= J2i\u a (r)\ 2 5(E - E a ) + \v a (r)\ 2 8(E + E a )}. (26) 



In order to calculate various physical quantities, we 
must construct the Green functions from E a , u a (r), 
v a (r) defined as 



)(x,x ) 



t\ ( 9n(x,x') gi2{x,x') 



g2i(x,x') g 2 2(x,x') 



- (T T [<h(x)ty(a/)]) - (Tri^ix^ix')] ' 



(18) 



- (TS\{x)^\{x')]) - (Tr^lix^ix')} 

with x = (x, t) . After the Fourier transformation of r as 



For finite temperatures, the (5-functions in eq. ( p6| ) are 
replaced by the derivative f'(E) of the Fermi distribution 
function f{E): 

N(E, r) 

= -Y}W a {v)\ 2 f{E a -E) + \v a (r)\ 2 f(E a + E)). (27) 



This finite temperature LDOS corresponds to the differ- 
ential tunnel conductance of STM experiments. 
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E. Nuclear relaxation time 

We now evaluate the spin-spin correlation function 
X-,+ (r, r', zi7„).Eil Its derivation is given in Appendix 
We obtain the nuclear spin relaxation rate by using Eq. 

§§, 

R(r, r') = Imx-,+(r, r', iQ„ -» Q + i V )/(n/T)\^ Q 

= -J2 u a (r) U ;(r)[u a (r')u* (r') + ^(r')^(r')] 

a,/3 

xnTf'(E a )6(E a -E p ). (28) 

We consider the case r = r' by assuming that the nuclear 
relaxation event occurs locally such as in Cu-site of high 
T c cuprates. Then, the r-dependent relaxation time is 
given by 

Ti(r) = l/JJ(r,r). (29) 

In Eq. d28|), we use S(x) = 7T _1 Im(;r — ir/) _1 to consider 
the discrete energy level of the finite size calculation. We 
typically use -q — O.Olt. 



F. Parameters used and site assignment 

The following parameter values are chosen: The aver- 
age electron density per site ~ 0.9 by adjusting appro- 
priately the chemical potential value \i. We normalize 
all the energy scales by the transfer integral t. For the 
s-wave case U = — 2.32< and V = 0. The resulting order 
parameter Ao/t = 0.5 at T=0 and H=0. For the <i-wave 
case U = and V = -A.2t. Then A /t = 1.0. The vor- 
tex unit cell is shown in Fig.[l], where the selective sites 
(V,A,B,C, and S-sites) to consider N(E, r) and Xi(r) are 
also indicated. The magnetic field H is denoted by the 
unit cell size N r as Hn t . Thus as N r increases, H itself 
decreases. 



III. BASIC VORTEX PROPERTIES 

In order to facilitate the understanding of the Mis- 
behaviors in the later section, we show some of the basic 
properties of vortex lattice. The following results coin- 
cide basically with those in the previous our own calcula- 
tions either in quasi-clas^icaiihepry for the vortex lattice 
of s- and d-w&ve cases Jl3Jiita~EZl and in BdG for isolated 
vortex of the s-wave casellj'O and-also with the preceding 
work by Wang and MacDonald.Ej 
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FIG. 1. A quarter region of the vortex unit cell which 
contains two vortices. Position of the sites V, A, B, C, and 
S in the square vortex lattice is indicated, where the nearest 
neighbor vortex is located in the 45° direction from the a-axis. 
The vortex center is shown by 0. The solid lines show the 
square atomic lattice. 

A. Order parameter profiles 



The pair potential or the order parameter A(r) van- 
ishes at the vortex center. |A(r)| exhibits the Friedel 
oscillations around the core both for the s-wave and d- 
wave cases whose period ~ {hp is the Fermi wave 

number). The amplitude of this quantum oscillation in- 
creases as the attractive interactions \U\ and \ V\ become 
large because the quantum effects are enhanced when 
Aq/Ef increases (E-p is the Fermi energy). These char- 
acteristics coincide with those in the s-wave case in our 
previous study.E2l 

As for the temperature dependence of the order pa- 
rameter, the vortex core radius, shrinks with decreasing 
T by the Kramer-Pesch effect MB We confirm it for the 
d-w&ve case too. The shrinkage is saturated at a low 
temperature both in the s- and d-wave cases. There, 
the structure of A s (r) and A<j(r) is almost independent 
of T . This is a-quantum-limit effect which occurs for 
T/T c < A Q /E F M We calculate the low temperature be- 
havior of T\ (r) by using the saturated pair potential. At 
higher temperatures, we calculate Ti(r) by using the self- 
consistently obtained pair potential at each T. 
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FIG. 2. (a) Local density of states N(r,E) for the s-wave case at T — and -H32 along the NN direction (V — > S). The 
backside (front) corresponds to V (S)-site defined in Fig. [lj (b) Local density of states N(r, E) along the NNN direction (V — * 
C). Note that the foot at the lowest levels at E/Aq ~ ±0.28 extends outward (arrows) along the NN direction. 



B. Local density of states around a core 

The LDOS is displayed as a function of the spatial po- 
sition for the s-wave case in Fig. || and for the d-wave 
case in Fig. || at H^- The spectral evolutions in Fig. 
U are shown for two different paths; along the nearest 
neighbor (NN) direction [Fig. ^(a)] and the next nearest 
neighbor (NNN) direction [Fig. [§(b)]. The low-lying ex- 
citations at the core, which correspond to the prominent 
peak structure, are characterized by the angular momen- 
tum (see for details Ref. |l^). The higher and higher an- 
gular momentum states situated at the higher energy are 
activated successively when going outward, forming the 




E/Ao 

FIG. 3. Local density of states N(r, E) at T = and H 32 
along the NNN direction (V — ► C) for the d-wave case. The 
backside (front) corresponds to V (C)-site defined in Fig. |lj 
The spectral evolution along the NN direction, which is not 
shown, shows similar behavior to this figure. 



characteristic spectral evolution in Fig. g. The lowest 
excited state which is seen as the highest peak extends 
towards the NN direction rather than the NNN direction 
as indicated by arrows in Fig. ||(a). 

The LDOS for the d-wave case is shown in Fig. |[ Since 
there is no discretized bound state around the core in this 
case, only the broad resonance peak centered at E = 
is seen at the core, which evolves smoothly. At the far- 
thest site the spectrum exhibits the known bulk behavior 
expected in the d-wave superconductor, characterized by 
_E-linear behavior due to the line node. In contrast with 
the s-wave case, there is no clear directional dependence 
in the spectral evolution in this scale (also see Fig. 0). 




E / A o 

FIG. 4. Local density of states jV(r, E) at T = and H Z2 
along the NNN direction (V — * C) for the d-wave case. The 
backside (front) corresponds to V (C)-site defined in Fig. [l]. 
The spectral evolution along the NN direction, which is not 
shown, shows similar behavior to this figure. 
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FIG. 5. Local density of states N(r, E) at T = and H32 along the NNN direction (V — » C) for the d-wave case. The 
backside (front) corresponds to V (C)-site defined in Fig. [j]. The spectral evolution along the NN direction, which is not shown, 
shows similar behavior to this figure. 



The spectral evolution in a high field Hn is depicted 
for the d-wave case in Fig. ||. The quantum oscillations 
due to the Landau band quantization becomes evident at 
a high field. We see rapid oscillations outside the main 
gap (\E/Aq\ > 1.0) and slow oscillations inside the main 
gap (\E/Aq\ < 1.0). Note that the split features around 
E ~ at the V-site in Fig. [| is not due to the localized 
bound state formation. In order to further facilitate the 
understanding of the T-dependence of the site-dependent 
Ti(r) behavior, we look into the spectral evolutions for 
some more details at H20 where we calculate T\ (r) in the 
next section. The LDOS around the vortex for the s-wave 
case is shown in Fig. ^(a), and the d-wave case is shown 
in Fig. |B|(b). In N(E, r) at the vortex center (the V-site), 
the gap edge at A in the zero-field case (dotted line U) 
is smeared, and low-energy peaks of the vortex core state 
appear. In the s-wave case, we see some peaks above the 
small gap Ai (~ Aq/Ep). It is due to the quantization 
of the energy level in the s-wave case. In the d-wave 
case, the core state shows zero-energy peak instead of 
the split peaks in the s-wave case.li3 There is no small 
gap. The weight of the low-energy states is decreased 
with going away from the vortex center (V^A^B^C). 
Far from the vortex, N(E, r) is reduced to the DOS of 
the zero-field case. But, small weight of the low-energy 
state extending from the vortex core remains there. It is 
noted that the weight of the low-energy state at the S-site 
is larger than that of the B-site in the s-wave case, while 
the S-site is farther from the vortex center [see lines for 
the S- and B-sites in Fig 0(a)] . It is due to the vortex lat- 
tice effect. The quasi-particle transfer between vortices 
occurs along the line connecting NN vortices (i.e., near 



the S-site) . This is also seen from Fig. ||(a) indicated by 
arrows. 

To obtain another perspectives of the LDOS image in the 
d-wave case we plot the spatial distribution of N(E = 
0, r) at H 32 in Fig. [| which directly relates to the 
Ti(r) behavior as discussed later. Notice that every- 
where N(E = 0,r) ^ 0. The LDOS N(E = 0,r) is 
largest at the core site and extends towards the NN di- 
rection. This shows the characteristic quasi-particle tra- 
jectory in the d wave case. The NMR experiment can im- 
age this low-energy excitation as is explained later. The 
detailed expnaitions of thc_quasi-particle trajectories in 
the s-waveoEj and d-wavg^3 are given previously. 




FIG. 6. Topographic view of the local density of states 
at the Fermi level normalized by the normal state value: 
N(E = 0,r)/jV(0) for the d-wave case at -H32. One unit cell 
(32x32 atomic sites) is shown. The vortices are located at 
the center and corners of the figure. 
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FIG. 8. Histograms p(H) of the amplitudes of the magnetic field for the s-wave case (a) and the d-wave case (b). 
should correspond to the NMR line shape. The letters in the figures denote the sites defined in Fig. |l| 



This 



C. Field distribution in vortex state 

In order to determine the one-to-one correspondence 
between the site and the magnetic field, we calculate the 
field distributions for the s-wave and d-wave cases. The 
magnetic field is evaluated via the current distribution 
through the Maxwell equation: 



V x H(r) = — j(r), 

c 



(30) 



where the current j(r) is calculated as 

jt{Ti) = 2\e\dm{t i+Sii J2{H(v i+i )^ a (v i )}} (31) 

= 2\e\clm{t i+ ea^2[u* a (r i+ e)u a (ri)f(E a ) 

(X 

+v a (r i+s )v* a (r i )(l-f(E a ))]} (32) 



for the e-direction bond (e = ±i, ±y) at the site r^. 

Figures. 0(a) and (b) show the contour plots of the 
field distributions for one unit cell at H^- It is seen 
that for the s-wave case the contour plot exhibits rather 
square-like shape around the core, which is contrasted 
with a cylindrical shape in the d-wave case. Due to the 
inter-vortex interaction, the magnetic field around a vor- 
tex tends to extend to the NN vortex direction. Then, 
the field of the s-wave case is modified to a square shape, 
reflecting the square vortex lattice shape. In the d x 2_ y 2- 
wave case, the magnetic field tends to extend tOpS-paida 
and b axis directions due to the gap node effect 
Then, the inter- vortex interaction effect and the gap node 
effect cancel each other, resulting in a cylindrical shape. 

The internal field distributions as a function of the 
magnetic field, which correspond to the resonance fre- 
quency distribution in NMR experiments, are depicted 
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in Fig. ||. We identify each site (A, B, C, S, and V) in 
the vortex lattice to these distributions. The one-to-one 
correspondence between the site position and resonance 
frequency can be used as a guide for site-selective NMR 
experiment: Namely, as shown in Figs. || (a) and (b), 
the NMR signal at the maximum cutoff of the resonance 
spectrum as a function of applied field or probe frequency 
comes from the vortex core at the V-site. With going 
away from the center (V— >A— >B— >C), the resonance field 
is decreased. The signal at the minimum cutoff comes 
from the C-site. The logarithmic singularity of the res- 
onance field comes from the saddle point of the field at 
the S-site. Thus it is possible to perform the site-selective 
Ti (r) measurement by tuning the resonance frequency. 

IV. TEMPERATURE DEPENDENCE OF Ti 

A. s-wave case 

We now consider the T-dependence of Ti(r) at each 
site defined in Fig. [I], which reflects the LDOS discussed 
above. The s-wave case is shown in Fig. ^|. We plot 
Ti(r) -1 vs. T for each site in Fig. ||(a), and re-plot it 
as lnTi(r) vs. T -1 in Fig. ||(b). We also calculate the 
zero-field case in our formulation. It is shown as line U 
in the figures. At the zero field, Ti ~ e A °/ T . Then, the 
slope of the InTi vs. T -1 plot gives the superconducting 
gap Ao, as the line U in Fig. ||(b). In the presence of 
vortices, Ti deviates from the relation e A °/ T at low T 
due to the low-energy excitation around the vortex core. 
This deviation was reported in the cxperiments.u In our 
results, reflecting the small gap Ai in the s-wave case 
[see Figs. ||(a) and (b), also see the V-site in Figs. 0(a)], 
Ti shows the slope Ai at low T in the InTi vs. T plot 
as seen in Fig. |(b). That is, T x ~ e Al / T . With leaving 
the vortex center, since the amplitude of the low-energy 
bound states is damped, the weight of e Al / T gradually 
decreases. Then the crossover temperature from e A °/ T 
to e Al / T is lowered. It is noted that Ti is faster at the 
S-site than that of the B-site, while the S-site is further 
from the vortex center. This non-trivial result is due to 
the vortex lattice effect noted above [see the foot denoted 
by arrows in Fig. ||(a)] . We should also notice the behav- 
ior of the coherence peak below T c . As seen in Fig. ||(a), 



with approaching the vortex center as C— >B, the coher- 
ence peak is suppressed. But in the vortex core region 
(lines V and A), a large new peak grows at intermediate 
temperatures. This is because the LDOS at the vortex 
core has peaks at low energy Ai instead of the singularity 
of DOS at Ao as seen from Figs. ||(a) and (b). 

B. d-wave case 

As for the d-wave case, we plot Ti(r) -1 vs. T in Fig. 
|Io|(a), and re-plot it as a log-log plot in Fig. |To|(b). 
At zero field (line U), we see the power law relation 
Tf 1 - T 3 of the d-wave case as expected. Note that 
this can be seen only below T/T c ~ 0.1 in our case. In 
the presence of vortices, Ti(r) -1 deviates from the T 3 - 
relation, and follows Ti(r)~ 1 ~ T at low temperatures. 
This deviation was reported in the experiments on high- 
T c cuprates.Q The origin of the T-linear behavior is the 
low-energy state around vortices in our case, instead of 
the residual density of states due to impurities or de- 
fects. With approaching the vortex center, the T region 
of the T-linear behavior is enlarged and it appears from 
higher temperatures. As seen in Fig. ||(b) of the d-wave 
case, the superconducting gap is buried by the low-energy 
state around vortices without the small gap of the order 
Aq/Ef- Then, T^ 1 ~ T at low temperatures in the d- 
wave case instead of the relation Ti ~ e Al / T in the s-wave 
case. As seen in Fig. [lO], Ti(r) -1 at the vortex center 
(line V) is very large compared with the zero-field case 
(line U). It reflects the fact that the LDOS of the low- 
energy state is larger than the DOS of the zero-field case 
as seen in Fig. ||(b). Figure |^ also shows the LDOS at 
E ~ normalized by N(0): the DOS in the normal state 
at the Fermi surface. It clearly indicats that the LDOS 
around the core exceeds the normal state value. This 
short relaxation may be the evidence of the low-energy 
peak in the LDOS by the low-energy core state. 

The coherence peak below T c is taken as a manifes- 
tation of the s-wave symmetry. In the d-wave case, the 
coherence peak is absent. But in the vortex core region, 
T^ 1 has a peak below T c even in the d-wave case. We 
should be careful not to mistake this peak due to the 
vortex core relaxation as the usual coherence peak in the 
NMR experiment when identifying the gap symmetry. 
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FIG. 9. Temperature dependence of Ti(r) in the s-wave case at the sites V, A, B, C, and S assigned in Fig. |lj. (a) 
Ti(T c )/Ti(T) is plotted as a function of T/T c . (b) lnTi(T)/lnTi(T c ) is plotted as a function of T c /T. Line U shows the zero 
field case. The line N is for the normal state at T > T c . 




V. FIELD DEPENDENCE OF T x 

So far, we only consider Ti(T) under a fixed H{— H20), 
focusing on the T-dependences both in the s-wave and 
d-wave cases. In this section we investigate the field de- 
pendence of the Ti(T) behavior. 



A. Field dependence of Ti(T) 

We calculate the T-dependence of Ti (T) for three rep- 
resentative fields; the weak field H32, the intermediate 
field H20, and the high field f/12 in order to see how 
the coherence peak evolves and the overall T-dependcncc 
varies with field. We show the results for the V-site and 
the S-site (other sites exhibit a similar behavior to the 



S-site, thus not shown here). They reflect the field de- 
pendenca-of the LDOS at the vortex core and outside of 
the core.Ej 

First, we consider the s-wave case. It is shown in Fig. 
|l2| . As seen in Fig. [l2](a), the T^ 1 - enhancement at the V- 
site in the intermediate temperature diminishes quickly 
with increasing field. This is because the sharp low en- 
ergy peaks of N(E, r) at the vortex core, which is the 
origin of the enhancement, is smeared by the effect of the 
quasi-particle transfer between vortices. As shown in Fig. 
|l2|(b) for the S-site, the coherence peak below T c also di- 
minishes with increasing field. In particular, for the high 
field H12, T^ 1 shows a depression below T c rather than 
the enhancement. This is caused by the smearing of the 
superconducting gap by the low energy state around the 
vortex, which extends outside of the core region by the 
quasi-particle transfer between vortices. With increas- 
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ing field, the singularity of the LDOS at the bulk gap 
edge E=Ao is smeared. Then the coherence peak dimin- 
ishes. Jt is qualitatively consistent to the observation on 
V 3 SnS 

The similar field-dependence of T\ occurs in the d-wave 
case. The T-dependence of Tf 1 is shown in Fig. [l3| At 
the V-site [Fig. |l3|(a)], the enhancement of Tf 1 is de- 
pressed as H increases. It reflects that the low energy 
peak of the LDOS around the vortex is smeared by the 
quasi-particle transfer between vortices. At the S-site 
[Fig. l^(b)], Tf 1 below T c is suppressed as H increases. 
It reflects that the low energy state around the vortex 
smears the d-wave superconducting gap-edge. The log- 



log plots of 7\ 



T are presented in Fig. U3 to see 



the low temperature behavior. At the V-site [Fig. |l3|(a)], 
the T-linear coefficient in T-f" (T) is depressed with in- 
creasing H at low temperature. This is because the 
inter- vortex quasi-particle transfer smears the low energy 



quasi-particle peak of the LDOS at the vortex core. At 
the S-site [Fig. ||(b)], Tf \T) is increased with raising 
H at low temperatures, while Tf ^T) is decreased with 
raising H near T c [Fig. |l^(b)]. Because the vortex contri- 
bution is increased and the amplitude of the low energy 
state extending outside the vortex core becomes large, 
the relaxation time at low temperatures becomes short 
with increasing external magnetic field. This tendency 
coincides qualitatively with the observation of high-T c 
cuprates by Ishida et alEs or an organic superconductor 
K-(ET) 2 Cu[N(CN) 2 ]Br by Mayaffre et alH 

We note that the quantum oscillations due to the Lan- 
dau band quantization affects the low temperature be- 
havior of Ti at extreme high fields. In Fig. (y| Tf 1 (T) 
at the high field H12 deviates from the T-linear behavior. 
It deviates downward (upward), when the LDOS is min- 
imum (maximum) at the Fermi energy in the quantum 
oscillation of the LDOS (Fig. |). 
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FIG. 11. Ti(To)/Ti(T) is plotted as a function of T/T c for the s-wave case at H32, H 20 , and H 12 . T c and Ti(T c ) are defined 
by each H. (a) the V-site and (b) the S-site. 




FIG. 12. Ti(T c )/Ti(T) is plotted as a function of T/T c for the s-wave case at H32, #20, and H 12 . T c and 2i(T c ) are defined 
by each H. (a) the V-site and (b) the S-site. 
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FIG. 13. The low temperature behavior of Ti(T)" 1 . The 
d-wave case at H32, H28, H24, H20 and the high field J/12. T c 
S-site. 

B. 7i(r,T) vs. LDOS 



Here we investigate the relationship between T 1 _1 (r) at 
low temperatures and the low energy excitations N(E = 
0, r) in the d-wave case. In the low temperature region, 
Ti shows T-linear behavior. First, we plot the spatial 
distribution of T 1 _1 (r,T ~ 0) at if 32 for a unit cell in 
Fig. [TJ. This shows that Tf l (r, T ~ 0) normalized by 
the normal state value T± is largest at the core and ex- 
ceeds its normal state value. Namely, with approaching 
the vortex core, T£~ increases. This tendency is in qual- 
itative agreement with a measurement by Milling and 
SlichterEl on 63 Cu NMR experiment of YBCO. 

The histogram p(Tf x ) of T x _1 (r,T ~ 0) is shown in 
Fig. where the values of Tf 1 (r,T ~ 0) for the to- 
tal 32x32 sites is classified according to its magnitudes. 
This indicates that the vast majority sites outside the 
vortex core exhibit long relaxation times, but certain few 
sites around the vortex core exhibit short relaxation time, 
which contains useful information on the site-dependent 
low energy excitations associated with the vortex. 

This topography of T 1 _1 (r, T — 0) in Fig. Wj^ looks sim- 
ilar to that of LDOS N(E = 0,r) in Fig. § In fact, as 
shown in Fig. |TJ where T^{t,T - 0) vs. N 2 (E = 
0,r)/N 2 (0) is plotted, the linear relationship between 
them is apparently seen. It implies that T 1 ~ 1 (r, T ~ 0) at 
low T is governed by the low energy excitations at each 
site. Then, the site-selective T 1 _1 (r, T ~ 0) is a good 
measure of the local core excitations. The LDOS around 
the vortex core can be estimated quantitatively from the 
slope of Tf^r) at low T. 



log-log plots of Ti(T e )/Ti(T) vs. T/T c are presented for the 
and Ti (T c ) are defined by each H. (a) the V-site and (b) the 
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FIG. 15. Histogram p(T 1 ~ 1 ) of Tf 1 corresponding to Fig. 

E 
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{N(E~0,r)/N(0)} 

FIG. 16. The linear relationship between T? /Ti(r) and 
[N(E = 0,r)/iV(0)] 2 . We plot the 32 x 32 latticepoint val- 
ues of the normalized T 1 v /Ti(r) shown in Fig. |lj and the 
square of the normalized LDOS N(E = 0, r)/7V(0) shown 
in Fig. |(| The letters in the figures denote the sites de- 
fined in Fig. [jj The dotted line shows the linear relation 
T 1 N /T 1 (r)oc[N(E = 0,r)/N(Q)] 2 . 



VI. CONCLUSION AND DISCUSSIONS 

We have calculated the nuclear relaxation time Ti in 
mixed state in the microscopic framework of the BdG 
theory and demonstrated that the site-selective Ti is a 
good probe to extract detailed information on low ly- 
ing vortex core excitations. We have investigated the 
site-dependence, T-dependence and field-dependence of 
Ti both for the s-wave and d-wave superconductors. 

Traditionally, the vortex contribution was considered 
as tbf|-spin diffusion to the normal region of the vortex 
core,Ht3 and Ti is treated as the spatial average. How- 
ever, we can investigate the position dependence of Ti(r) 
around vortices through the resonance field dependence. 
This is an advantage of NMR over other methods. We 
should clarify the local mechanism of the relaxation (i.e., 
whether the relaxation occurs locally, or it is averaged 
by the spin diffusion). It is noted that in the clean limit 
the vortex core regitm-is- noj . a simple core filled by nor- 
mal state electrons J'EraOETEZl There, the characteris- 
tic T-dependence is expected near the vortex core other 
than a simple T-linear behavior, reflecting the rich struc- 
ture of the low energy state around the vortex core. 

We expect that the NMR imaging study just ex- 
plained here will provide vital information for the vor- 
tex core state in high-T c cuprates. As for the problem 
whether the quantization of the energy levels occurs or 
not, Ti ~ e Al/T if the gap Ai (~ Al/E F ) is present in 
the excitation due to the quantization. If this small gap 
is absent, T-j - ~ T. As for the problem whether the 
zero-energy peak exists or not in the core state, the re- 
laxation at the core becomes eminently faster than that 
of the zero-field case (or that far from the vortex) at 
low temperature, if the zero-energy peak exists in the 



LDOS as suggested in the theoretical study. If the peak 
structure is absent within A as reported in the STM 
experiments on BSCCO, the relaxation is slow even at 
the vortex core as in the zero-field superconducting case. 

We have pointed out a pitfall of the conventional pro- 
cedure of the NMR experiment to extract the nodal 
gap structure by analyzing the power law of the T- 
dependence of Ti or the coherence peak. Even in the 
d-wave pairing, Tj -1 shows enhancement below T c by the 
vortex core contribution. 

Finally, we stress that, to perform an idealized site- 
selective NMR measurement, we need a clear NMR reso- 
nance line shape of the vortex lattice as shown in Fig. 
^ or as that obtained by a beautiful experiment on 
vanadiumfl'D In the resonance line shape, the signal near 
the upper (lower) cutoff of the field distribution gives 
the information from the vortex core (the outside of the 
core). However, even in the case when the resonance line 
shape of the vortex lattice is not clear, if we analyze the 
fast and slow relaxation processes separately, the fast re- 
laxation process includes the information of the vortex 
core contribution.EJ We hope that future NMR experi- 
ments confirm our proposal of the site-selective NMR, 
and clarify the relation of Ti(r) and N(E,r). If this 
experimental method is established, it can be a powerful 
method to investigate the exotic mechanism of the uncon- 
ventional superconductors by spatially imaging the low 
energy quasi-particle excitation around the vortex cores. 
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APPENDIX A: 

The equation of motion for the Green function in a 2 x 2 matrix form is written in the coordinate space as 

£{-ft|-l*(x - x") - ( f/*' _^ X X 2 ) }3(x"r, xV) = M(x - x')5(r - r')l (Al) 

x" 

where 1 is the 2x2 unit matrix. We put the completeness relation: 

u a (x) -<(x) \ / <(x') <(x') 



E a >0 



V 7 ^ 1 D„W <( X ) / V ~^( x ) u a( x ) / 



before j(x"t,x't') in Eq. (Al), and use the BdG equation of the form 



v / if(x,x') A(x,x') \ / u a (x') -<(x') \ _ / u a (x) -<(x) \ / E a \ 



(A4) 



1> V At ( x ' x ') ~K*{^) ) \v a (x.') u*(x') J \v a (x) <(x) ^0 -£ a 
for E a > 0. After the Fourier transformation of t to u„, the inversion of the matrix product leads to 

<Kx,x',o,„)= y ( Ua{ *\ - V °W)( ? ) ( <{ f\ 



with the help of the completeness and the orthogonal relation of the wave functions. In the term with (iuj n + E a /K) 1 , 



we use the symmetry property: u a — > — w*, v a — > it* when E a — > — E a . Then, Eqs. (20)-(E3J) are obtained 



APPENDIX B: 

The spin-spin correlation function is defined as 

X-.+foz 7 ) - (T T [S-(x)S+{x')]) = (T^lix^ix^ix'^^x')}) (Bl) 

with x = (x, t). This can be rewritten in terms of the thermal Green functions, 

X-,+ (x,x') = g 11 (x,x')g 2 2(x,x') - g 12 (x, x')g 2 i(x, x'). (B2) 
After the Fourier transformation of r, we obtain 

X-,+ (x,x / ,O n ) = T^{gii(x,x',w„)522(x,x',fi„ - u n ) - ,gi 2 (x, x', w„).g 2 i(x, x', ft n - w n )}. (B3) 

Then, substituting Eqs. (p0|)-(p3|) and performing w n -sum, the spin-spin correlation function with real frequency Q is 
reduced to 

X -, + (x,x',iQ n -> n + lV ) =-J2 u»(x)^(x)«(xO^(xO + v* a (x')v fj ( X <)} f{Ea) - f{E ?\ . (B4) 
Here, we used the symmetry property: u a — > — v* , Uq, — > it* when _E Q — > —E a . 
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